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Abstract 

By using the Jacobi metric of the configuration space, and assuming ergodicity, we cal- 
culate the Boltzmann entropy S of a finite-dimensional system around a non-degenerate 
critical point of its potential energy V. We compare S with the entropy of a quantum 
or thermal system with effective potential V e ff. We examine conditions, up to first 
order in perturbation theory, under which these entropies are equal. 
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The Boltzmann entropy S is a thermodynamic potential which is particularly useful 
in describing the macroscopic behavior of a system of constant energy E. We calculate 
the entropy variations due to quantum or thermal effects of a finite dimensional Hamil- 
tonian system 1 ' 2 , where we assume that these fluctuations are encoded in the effective 
potential 3 ' 4 ' 5 V e ff of the system. S is expressed as the logarithm of the volume of the 
configuration space M. of the system under study. Calculating such a volume is im- 
practical for most systems, so to simplify the problem we confine ourselves to spheres 
around a critical point P of the classical potential V c . The evolution of a physical 
system in equilibrium takes place in a neighborhood of a local minimum of V c compat- 
ible with the external conditions, hence the importance of computations taking place 
around such a critical point. 

Determining volumes on a manifold does not require the introduction of a metric. 
The existence of a top- dimensional differential form is sufficient, as in the case of the 
Liouville form on phase space. However, if the manifold possesses a metric, especially 
a Levi-Civita one, considerable simplifications occur due to the strong constraints im- 
posed upon it. On the other hand, a physical model depends to a great extent on the 
metric chosen for analyzing a system. Many times there is no a priori justification of 
the choice of a specific metric. For our analysis we follow essentially Krylov's ideas 6 
and use the Jacobi metric g, which has components in a coordinate basis 

g ij = 2[E-V(x)]5 ij (1) 

where E is the total energy of the system and V(x) is its potential. From a mechan- 
ical viewpoint, the choice of this metric is convenient since it is directly related to 
the kinetic energy of the system. Another advantage of the Jacobi metric is that it 
is reasonably easy to manipulate, since it is diagonal and explicitly conformally flat. 
Because of its relation to the total energy of the system E, g is ideal for analyzing 
a system described macroscopically through its entropy S. The only disturbing fact 
about the choice of g is that, it is identically zero at the "turning points", namely at 
subspaces of the configuration space at which the kinetic energy is zero. However, if 
all calculations take place away from the "turning points" , g can be used reliably to 
analyze the behavior of the system. In the present paper, using g we find that the con- 
dition for S to remain invariant under a perturbation, is that the perturbing potential 
should satisfy a "massive" elliptic differential equation, and comment on its solutions 
in two extreme cases, at which considerable simplifications occur. 
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To begin with, an explicit expression for the Ricci scalar R of g is needed. The 
Christoffel symbols are given for a general metric, in this coordinate basis, by 7 , 



T )k = \d im {dj9mk + d k g jm - d m g jk ) (2) 



For g they become 



1% = - 2{ e_ v) + §W - 5 im 5 jk d m V) (3) 



The Riemann curvature tensor for any Levi-Civita metric is given, in this coordinate 
basis, by 

D 
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which, for g gives 
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(SlSidjVdiV - S^djV^V + 5}d k Vd m V + 

5%d m Vd 3 V - 5 a 6 im d k Vd l V - Si^SjmdiVdrV + ^Sj^diVdrV - 
S^dkVdiV + 6% m d k Vd t V - 5\d 3 Vd m V - 5%d m Vd l V + 
8 u 5 km d j Vd l V + 8)8 rl b km diVd r V - 5 m 5)5^8^ d n V) (5) 

Contraction of % and k gives the Ricci tensor components 

R " = ^If^l^i^/^l+^^Ml^^HM^Jvyi 2 ] (6) 

In (6) A, V and || • || are the Euclidean Laplacian, gradient and norm respectively. 
Here n is the dimension of the configuration space M.. We raise one index by using 
the inverse of the Jacobi metric g _1 , and upon one further contraction we finally get 
the Ricci scalar 

«=£^-<!g^||W||> (T) 

One immediate observation is that the Jacobi metric g is not an Einstein metric, since 
it is, obviously, impossible to find a constant A such that Rij = Kg i3 even at the critical 
points of V. This does not allow us to take advantage of the rather simple variational 
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characterization 8 of the Einstein metrics. 



Our computations take place around a critical point P of the classical/non-thermal 
potential V c (x). For simplicity, we also assume that P is non-degenerate. If it is 
degenerate, we can use equivariant extensions of these arguments, especially if the 
degenerate critical points are connected by a transitive, free group action. For P, 
equations (6) and (7) for the Ricci tensor and the Ricci scalar, respectively, simplify to 



and 



^ = ^E^V) [{n ~ 2)didjV+6ijA2V] (8) 
n — 1 

R = -T rAV (9) 

2(E-V) V ; 

We also present, for completeness, the calculation of the volume of the image of a ball 

of radius r in the configuration space 9 . Consider a normal coordinate system whose 

origin is at P. Let O p denote the origin of a corresponding coordinate system in T P M. 

We restrict ourselves to balls B n (O p ,r) of radius r in T P M for which the exponential 

map exp p : B n (O p ,r) — > B n (P,r) is a diffeomorphism. In other words, we restrict 

ourselves inside the injectivity radius of P. However, we have already assumed in our 

treatment that we want to be away from the "turning points" at which the metric 

vanishes degenerately. So, for practical purposes in our treatment, r is less than the 

minimum of the injectivity radius of P and the distance to the closest "turning point" 

from P. We parametrize the points of i 6 B n (O p ,r) \ {O p } of T p M in a polar 

coordinate system as 

(\\x\\,-f-)e(0,r)xS n ^ 

\\X\\ 

where S^-i is the unit sphere in n dimensions. Then the calculation reduces to 
determining the "radial" volume element in At 10 . Consider an ortho normal basis 
{e±, e2, . . . , e n _i, X}, where X is the unit vector, with respect to the canonical metric 
in (0,r), which essentially defines a point on S n -±. Let us consider the geodesic c(t) 
on At with tangent vector X and the transversal Jacobi fields Ji(t), ^(i), . . . , J n -i(t) 
that satisfy Jj(0) = 0, -DJj(O) = e^, % — 1, . . . , n — 1. In this notation D stands for 
the unique Levi-Civita connection on TM. compatible with g. The "radial" volume 
/x t (X) is given by the absolute value of the Jacobian of Jj, i — 1, . . . n — 1 or since Ai 
possesses a Riemannian metric by 

d[j, t (X) = y/det g(J h J k ) dt, i, k = 1, . . . , n - 1 (10) 



3 



Moreover Jj(t) = tD[exp(tX)ei\ and since {e^}, i — 1 . . . , n — 1 is an orthonormal 
basis, we find that 

dfit(X) = * n_1 v / det^(exp(tX),exp(tX)) (11) 

Obviously = * n_1 dt for i = 0, which expresses the fact that the Riemannian 

manifold, in normal coordinates, is locally Euclidean. The volume of B n (P, r) is given 

volB n (P,r)= [ dS n -! f dfi t (X) (12) 

JSn-l JO 

In normal coordinates, it is well-known 7 that around point P G M. the metric can be 
expressed in terms of the Riemann curvature tensor and its covariant derivatives as a 
power series with lowest order terms given by 

g tJ (x) = % + l -R ikjm {P)x k x m + 0{\\xf) (13) 

We recall 11 , that the sectional curvature K along a plane spanned by the linearly 
independent vectors Y = Y % Ui and Z = Z % u,i is a function on the Grassmann manifold 
G 2 ,n{M) given by 

{gikgjm - 9 lJ g k m)Y i ZJY k Z m 

Here i — 1, . . . ,n is a basis of T p Ai, not necessarily orthonormal. The Ricci 

curvature in the direction of Y — Y l Ui is the average of the sectional curvatures along 
the planes containing Y, namely 



Using (13), we find 



det 9ij (tX) = 1 - l - P {X)t 2 + 0(t 3 ) (16) 



so the infinitesimal radial volume element dfj, t (X), which is proportional to ^det 
becomes in this approximation 

diH(x) = r- 1 - \ P {x)t n+1 + o(t n+2 ) (17) 

o 

Therefore, the volume of the ball B n (P, r) of M. is 

volB n (P,r)= [ dS n _i I (t n - 1 - l -p{X)t n+1 )dt (18) 
Js n -i Jo b 
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This non-trivial integral essentially calculates the "average" of the Ricci curvature over 
all directions. We expect it, therefore, to be proportional to the Ricci scalar R at P. 
To formally justify this, we consider an orthonormal basis {e[, . . . e^^} in which 
the Ricci curvature, considered as a quadratic form, is diagonal. Then the elements 
e-, i — 1, . . . n — 1 of this basis are eigenvectors of the Ricci curvatures with eigenvalues 
Si, i — 1, . . . n — 1. In this basis, 

x = Y. xl < ( 19 ) 

i=l 

so the Ricci curvature and the Ricci scalar, respectively, are expressed as 

rt— 1 rt— 1 

p(X) = 5>pP) 2 , flp = J> (20) 
i=i i=i 

The volume u n of the unit ball in Euclidean space is given, in spherical coordinates, 
by 

w„=/ (X*) 2 ^-! (21) 

•'Er= 1 1 (* i ) 2 =i 

which gives 

p(X)dS n - 1 = R p u n (22) 



Moreover, 

r /•'■ " i 

dt (23) 



volB n (P,r)= [ dS n -! [ 

JS n -! JO 



t"- 1 - ^p(X)t n+1 



which becomes +2 

volB n (P, r) = -r n A n ^ - J R p u n (24) 
n 6(n + 2) 

It is also well-known 10 that the volume of a ball of radius r, is given for any n G N 
by 

n 

= TxfTT) (25) 

Since <iO„ = a n -\dr we find for the area of the unit sphere S n -\ that 

a n _i = nu; n (26) 

Putting (9), (24), (25), (26) together, we finally find for the volume of the geodesic ball 
B n (P,r) 

volBJP, r) = f{ — ) - 6( „ + 2)r( . + 1) *, (27) 
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The Boltzmann entropy is defined by 



S = ks lnvolB n (P, r) 



(28) 



where k B is Boltzmann's constant. This definition has physical significance as long as 
the system is ergodic, i.e. as long as the geodesic flow in M. of the system is ergodic. 
For such a system, combining (27) and (28) we find 



S = k B In 



T[2r 



n+2 



r(f + l) 6(n + 2)r(f + l) 



Let's assume that the renormalized effective potential of the system is 

V eff (x) = V c (x) + XV(x) 



(29) 



(30) 



Here A indicates a small coupling constant which is used as a parameter in a power 
series expansion. V(x) represents the quantum/thermal contribution to the classical 
potential. The corresponding Jacobi metric on M. is 

g ij = 2[E-V c (x)-\V(x)\5 ii (31) 
We observe that both metrics are conformally equivalent, and g can be written as 

V(x) 



gij = 2 [E - V c i 



x 



1- A 



E - V c (x) 



which means 



9ij 



1 - A- 



V(x) 



E - V c (x 

The conformal factor relating these metrics is therefore 



9ij 



(32) 



(33) 



a(x) = - In 



1 - A 



V{x) 



E - V c (x) 



(34) 



A laborious computation, done either straightforwardly in coordinates, or in a coordi- 
nate free way using Nomizu's approach 13 , gives 



where 



B] 



R = e -M*) [R-2(n-l)Bl] 



Q - l) gUfodta + g kl (d k d t a - T™d m a) 



(35) 
(36) 
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At the origin of the normal coordinate system, the Christoffel symbols are zero. There- 
fore B] at P becomes 



B] = 



, n 



A(E-V C -XV) 2 



[A 2 (- - Vg^VdtV - 2Xg kl (E -V c - \V)d k d l V 



V 



-2XVg kl {\ - X^y)d k d l V c + 2X 2 g kl d k Vd l V] 
Expanding in powers of A, and keeping up to linear terms in this parameter, we get 



B) 



X 



2(E-V C ] 



V 



g^dudtV + ——g^dkdtVc 
E — Vc. 



(37) 



This gives 



R p -R p = -2(n-l)- t 



X 



g kl d k d l v + —— g kl d k d l v t 



V 



kli 



2(E-V C ) 
We observe that to zeroth order in A 



E - Vc 



-X- 



V 



E-V 



Rp Rp 



R P (38) 



(39) 



so the difference between the entropies S and S of the original and the perturbed 
system, respectively, due to (29) becomes 



S-S = k R \n 



1 



Rp Rp 



6(n + 2) 1 



:R r , 



(40) 



6(n+2)-"-P. 

Substituting (38) into (40), expanding the logarithm in powers of A and keeping up to 
linear terms, we find 

„2 



S-S = k B (-X) 
n — 1 



6(n + 2)^ 



E — V, 



V 



6(n + 2) 



g d k d t V + — — —g d k diV c 

Hi — V r 



+ E^V Rp 



(41) 



In this approximation, the entropy remains invariant if 



n- 1 

E-V c 

which reduces to 



V 



g kl d k d { V + ——g^dkdtV 
& — v n 



+ 



V 



E-vr p 



R n = 



3 ^b^ + ^) A ^-° 



(42) 



(43) 
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which eventually gives 



AV(x) + ^jV(x) = 



(44) 



Let A g denote the scalar Laplacian of the configuration space with respect to the Jacobi 
metric g. Then 7 

A, = l-aibTy/gdj) (45) 
At P, we find straightforwardly, for V(x) 

AV(P) = 2(E - V c {x))A g V{P) (46) 
From (9) we can also see that the Ricci scalar at P can be expressed as 

R p = (n - l)A g V c (P) (47) 

This condition can be equivalently written in terms of the scalar Laplacian A g associ- 
ated with g, as 

3A,y c (P) 



A g V{x) + 



-V(x) = 



(48) 



2(E-V C (P)) 

Therefore invariance of the Boltzmann entropy under perturbations amounts to the 
perturbing potential V(x) obeying a "massive" Laplace equation with "mass" 

3\A g V c (P) 



M = 



2{E-V C {P)) 



(49) 



If P is a relative minimum of V c (x), then M is positive-definite. Indeed, P is a local 
minimum of V c (x) so the Hessian at P is 



HijV c (P) = didjV c (P) 
and is positive definite. By definition 8 ' 9,11 



trH i:j = A g 



(50) 



(51) 



so after diagonalizing Hij, we see that A g V c (P) > 0. 
A special case arises if 

E = + V C (P) 



n-2 



Then the equation for V(x) inside the geodesic ball becomes 



A + 



n-2 



R„ 



4(n-l) P 



V(x) = 



(52) 



(53) 
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which means that V(x) is an eigenfunction of the conformal LaplaciaU 
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(54) 



corresponding to zero eigenvalue. In case E has any other value, V(x) is a solution of 




This is very similar in form, though just a linearized equation, to the Euler-Lagrange 
equation arising in the variational determination of the Yamabe constant 15 of a con- 
formal class of metrics on a manifold M. The common point of these equations is that 
both arise from variations within a conformal class of metrics of A4. 

If E — > oo, then the Riemann tensor (5) and consequently the Ricci tensor (6) 
approach zero. In this limit M, is Euclidean and the induced Riemannian measure on 
it, i.e. the induced volume element, is uniform. Then the effective mass M in (48) is 
zero. This limit describes the evolution of a system of very high kinetic energy. Because 
of the high kinetic energy, as the system evolves in A4, it is unaffected by any curvature 
deviations of M. from flatness. For S to remain invariant under perturbations, to first 
order in A, the perturbing potential V(x) must be a harmonic function with respect to 
the Jacobi metric. Such a function depends on the boundary conditions imposed on 
the sphere dB n (P, r). If r is small enough, it is reasonable to assume that the value 
of V(x) at all points x of dB n (P,r) is the same, i.e. that 



The unique harmonic function satisfying this condition is the constant function V . 

The other limiting case is one in which E — > V C (P). Then the Riemann tensor (5) as 
well as the Ricci tensor (6) approach infinity, a fact that indicates the curvature tensors 
develop conical or cuspidal singularities, depending on the exact form of V C (P). The 
source of the singular behavior is that in this limit the Jacobi metric becomes exactly 
zero. Therefore it is not a reliable probe of the structure of M. and the description of 
the evolution of the system. In such a case, the system is evolving on M. around P 
having a very small kinetic energy. No curvature fluctuations of M. can be overcome by 
the system, exactly because of its small kinetic energy. Equivalently, we observe from 
(49) that M becomes very high. For the energy E to remain conserved the system has 



V(x)\ dBn (P, r ) = V 



(56) 
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to confine itself in a very small ball B n (P, r) of A4. In this limit, the Dirichlet problem 
implies that the fluctuating potential should vanish at dB n (P, r). The unique solution 
to the boundary-value problem is therefore 

V(x) = (57) 

so S remains invariant as long as the first order radiative corrections of V e ff are trivial. 

These two limiting cases reflect the fact that the only harmonic functions on a compact 
manifold are constants. For all other values of E, between zero and infinity the system 
interpolates between these two extreme cases. Depending on the value of r, R p and the 
sign of R p , the accessible part M. may be disconnected, or may even be topologically 
non-trivial. This fact places further restrictions on the form of V(x). 

When R p > there is no guarantee that the system is ergodic even on B n (P, r) 
for small r. Ergodicity of the geodesic flow associated to g is related to the Riemann 
curvature tensor (5) or equivalently to the sectional curvature K (14) of g rather than 
its Ricci scalar R p (9). There are various quantitative measures of the complexity of 
the geodesic flow of a manifold. One of the most studied measures is the topological 
entropy h t0 p 16 - It is known 16 that if the sectional curvature K (14) of a manifold M. is 
everywhere negative then its fundamental group has exponential growth 17 , the geodesic 
flow is ergodic and the topological entropy is positive. K is not negative, generically, 
for g on M. as can be seen from (14). Confining our attention to the sphere under con- 
sideration, it is impossible to have negative K everywhere inside B n (P, r) for any value 
of r. This is a direct consequence of Milnor's theorem 17 and the fact that iri(B n (P, r)) 
is trivial for n > 1. A metric of everywhere positive K however may or may not have 
htop > 0. Since K does not have a constant sign throughout B n (P,r), for the entropy 
S to be a physically relevant measure of disorder, we have to assume ergodicity of the 
geodesic flow of g. 

A by-product of these calculations is a contribution towards the computation of 
the Euler characteristic needed in the topological hypothesis 6 . According to it, the 
source of phase transitions lies in the topology change of appropriate sub-manifolds Af 
of M. This topology change is detected by calculating the Euler characteristic of Af 
usually by using Morse theory 6 . Let AT be a d— fold un-ramified covering of Af. Then 18 
x(AT) = dx(Af), so the calculation of the Euler characteristic of Af reduces to that of 



10 



AT where the integration measure in the Gauss-Bonnet-Chern formula 19 applied to J\F 
is provided by (27) as long as AT has the topological type of a homology m— sphere, 
m < n. 
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